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ABSTRACT 

Transmission  and  reception  of  electromagnetic  f  z^ ) 
energy  "by  communications,  weapons,  and  active/passive 
sensor  systems  is  known  to  be  strongly  influenced  by  an 
atmospheric  phenomena  known  as  ducting,  caused  tj 
refractive  layers  in  the  atmosphere  of  marine  environments. 
The  Naval  Postgraduate  School  (NFS'  has  developed  a  Marine 
Atmospheric  Boundary  Layer  (M.^3L)  model  which  can  r^  isc" 
to  predict,  over  a  24  hour  period,  the  refractivf  profile 
of  the  lower  atrrosphere.  This  thesis  examines  the  model 
from  the  statistical/trend  analysis  approach  to  examine 
whether  the  model  can  be  used  as  a  valid  predictor  of 
refractive/duct  i  n  ^  conditions  . 
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I.   INTRODUCTION 


The  phenomena  i n o w n  as  ducting  occur1  when  refractive 
layers  cause  electromagnet  ic  (EM)  energy  to  "rend  toward  the 
earth  at  a  rate  greater  than  or  equal  to  the  curvature  of 
the  earth.  Ducts  occur  in  two  forjis;  (1)  surface-based 
ducts  caused  by  either  a  surf ace-la sed  layer.  Figure  1,  or 
by  an  elevated  layer,  Figure  2,  and  (?)  elevated  du-ts 
caused  by  an  elevated  layer,  Figure  3.  As  shown  in 
Figures  1  through  3,  duct  types  and  thicknesses  are  easily 
determined  by  utilizing  the  profile  for  the  modified  index 
of  refraction  M  which  depends  on  temperature,  pressure,  and 
specific  humidity. 

The  Navy  employs  the  Integrated  Refractive  Fffects 
Prediction  System  (IEEPS)  to  identify  ducting  conditions 
and  to  assess  the  effects  on  various  fleet  emitters  [Ref. 
l].  A  drawback  to  I  REPS  usage,  from  a  tactical  standpoint, 
is  that  it  depicts  the  ducting  conditions  only  at  the  time 
of  the  radiosonde  sounding. 

From  a  tacticians  viewpoint,  the  ^-profiles  are  a 
valuable  output  of  the  Marine  Atmospheric  Boundary  Layer 
(MABL)  model  because  of  the  ease  with  which  duct  types  and 
thicknesses  are  obtained.  The  frequencies  which  are 
trapped   by   ducts   ar^  a  function  of   duct   thickness,   as 
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M  units 
Figure  1.  Surface  Duct  Caused  by  Surface-Based  Layer 
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Figure  2.  Surface  Duct  Caused  by  Elevated  Layer 
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M  units 
Figure  3.  Elevated  Duct  Caused  by  Elevated  Layer 


stated  by  Kerr  (1951).  Thereby,  knowing  the  duct  thickness 
and  type,  the  tactician  would  be  able  to  determine  the 
frequency  range  which  would  be  trapped.  In  order  for  the 
EM  energy  to  "be  trapped,  the  transmitting  antenna  must  be 
physically  located  within  the  duct.  Antenna  heights  are 
readily  obtainable  from  ship's  characteristic  cards,  hence 
the  transmitters  on  board  which  would  be  affected  by  the 
duct  could  be  quickly  identified.  The  ability  to 
accurately  predict  the  evolution  of  the  ducting  conditions 
for  a  24  hour  period  would  be  a  powerful  tool  in  the  hands 
of  today's  Carrier  rattle  Group  Commanders.  This 
information  could  be  used  to  establish  "mission  Control 
(EMCON)  policy  and  to  establish  effective,  efficient 
frequency  monitoring  plans.  The  MAEL  model,  developed  by 
the  Environmental  Physics  Group  at  the  Naval  Postgraduate 
School  (NFS),  may  be  the  medium  through  which  this 
information  can  be  made  available  to  the  fleet  tacticians. 

The  purpose  of  this  thesis  is  to  perform  a  detailed 
statistical/trend  analysis  of  the  MATL  to  determine  its 
validity  in  predicting  ducting  conditions  for  a  24  hour 
period.  Cnce  this  analysis  has  been  completed,  conclusions 
and  recommendations  will  be  made  in  order  to  provide  a 
foundation  for  further  refining  of  the  N'ABL  model,  if 
necessary,    by    the   Environmental   Physics   Croup.   The 


overall  objective  is  to  provide  a  reliable  predictive 
model,  available  for  dissemination  and  incorporation  in  the 
fleet  . 

The  approach  taken  in  this  study  differs  from  past 
MABL  model  evaluations  at  NFS,  in  that  this  approach  is  to 
perform  an  objective  examination  of  the  Todel  .  The  goal  is 
to  extract  information  about  the  model  output  from  a  purely 
statistical/trend  analysis  approach,  which  has  not  been 
d  one  previ  ously . 
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II.  ^Cr^L  INSCRIPTION 

T h s  v A B L  it o d e  1  is  a  ze:o-:rd?r  ,  two  layer,  integrated 
mixed  layer  no del.  The  atmosphere  Is  assumed  to  consist  of 
two  layers:  a  well-mixed,  turbulent  boundary  layer,  and  the 
relatively  non-turbulent  free  atmosphere  above.  The 
inversion  (or  "transition  zone")  separates  one  lever  frcrr 
the  other.  In  a  z e r o - o r d e  r  model,  t v e  i  n v e r s i  o r  thickness 
is  assumed  to  be  zero,  hence,  a  discontinuity  or  "jump" 
occurs  in  the  profiles  of  the  variables  at  the  inversion 
[Pef .  ?]. 

The  MPS  model  of  the  mixed-layer  is  ij^I-emente.d  on  a 
Fewlett  Packard  9845  desktop  computer.  The«mcdel  u^s  a 
30-minute  time  stej  to  predict  Tixed  layer  temperature, 
specific  humidity,  the  jump  of  these  values  at  the 
inversion,  and  cloud  or  fog  formation. 

The   M   profile  is  calculate!  from  predicted  T   and   c 

values  using  Equation  (1). 

5  2 

V   =  77.5(P/T)  +  6  x  10   *  (q  *  P/T  )  +  55.157  *  Z      (1) 

where  is  pressure  in  millibars,  T  is  temperature  in 
degrees  Kelvin.,  q  is  specific  rumidity,  and  Z  is  height  in 
meters.  The  top  of  the  durt  corresponds  to  the  height 
above   the  surface  where  the  M  value  is  minimum.    The  base 
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of  the  duct  is  the  heigh:  at  which  a  vertical  lice  irawn 
downward  from  the  point  where  1  is  a  minimum  value  to  where 
it  first  intersects  a  point  of  equal  ^-unirs  or  the 
surface,  whichever  sccurs  first.  Ducting  commonly  exists 
where  t here  is  a  temperature  inversion  which  acts  as  a 
trapping  layer,  refracting  F*l  energy  toward  the  earth. 
Inversions  typically  exist  in  marine  surface  high  pressure 
regions  where  well-mixed,  stable  layers  are  found  Vet  ween, 
warm  dry  air  above  and  cooler  moist  air  below.  Large  a-ras 
of  low  level  stratus  clouds  often  indicate  areas  of  iuct 
occurence.    Ducting  is  exuected  to  he  minimal  "ie  =  :   fronts 


and   areas  of  convective  cloudiness 


rents,   wit h   t  h e  i r 


associated  upward  motion,  often  dj<|4-pate  the  inversions  as 
the   whole  air  column  becomes  mixed,  f     Area1-:,  of   convective 

activity,   discernible   by  the  presence  of  cumulus   clouds, 
are  also  normally  inversion  free  [Pef.  1]  . 
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III.   CRITFFI_A  AND  ASSUMPTIONS  FOE  £ABL  ^OFEl  VEBI7ICATI0N 

A  major  assumption  in  the  verification  of  this  rod  el 
is  that  the  radiosonde  data  is  correct.  This  point  will  be 
elaborated  upon  later  in  this  paper. 

A.   CRITERIA  FOR  CHOOSING  EXAMINATION  PERIODS 

Forty  cases  (examination  periods^  were  chosen  to 
evaluate  the  MAEL  model.  The  radiosonde  soundings  used  in 
ail  4?  cases  were  taken  off  the  coast  of  Southern 
California  in  and  around  San  Clement?  Island.  Th^i^'- 
examination  periods,  which  arp  2*-  hours  in  length,  wf re 
chos pn  ba sed  on  the  following  criteria:  A 

1.  There  had  to  be  at  least  one  radiosonde  sounding 
in  the  26  hour  period  preceding  t'°e  examination  period  with 
surface  layer  data,  from  the  NFS  research  vessel,  Acania, 
within  2  hours  of  that  sounding.  This  requirement  is 
needed  to  calculate  Ws  (subsidence)  with  lenschcw's  (1973) 
method. 

2.  There  had  to  be  surface  layer  data  within  2  hours 
of  the  radiosonde  sounding  at  the  beginning  of  the 
examination  period,  also  required  in  calculating  Vs. 
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3.  There  had  to  re  at  least  2  soundings,  In  addition 
to  the  initial  sounding,  during  the  examination  period  to 
he  used  as  verifications  of  the  model's  predictions. 

4.  Wind  speeds  had  to  he  evailahie  within  3  hours  of 
the  initial  sounding.  If  wind  speeds  are  available  within 
3  hours  of  the  12  hour  point  or  within  3  hours  of  the  24 
hour  point  they  will  also  he  used.  If  these  2  values  J  re 
not  available,  the  initial  value  will  re  used  end  held 
constant  during  the  exam  period. 


B.   INITIALIZATION  REQUIREMENTS  FOR  M'OTEL  PREDICTIONS 

The  rode!  requires  as  input  data  the  following: 
1  •   E§dio sonde  S^und i rig  HsLla. 

Raw   data   from   the   radiosonde  launches  must  he 

digitized  for  use  in  the  VAPL  model.  ▼ 
2«   ilitsidence  iWs_l  ?redic_tion 

The  subsidence  value  used  in  the  model   prediction 


for   all   42   cases  was  obtained  from   integration 


the 


Moisture  Budget  Equation  given  in  Equation  (2).  This 
method  was  developed  and  performed  by  Lenschow  (1973).  It 
is  based  on  the  assumption  of  well-mixed  specific  humidity 
(q)  in  the  boundary  layer,  and  that  changes  occur  due  to 
fluxes  at  the  sea  surface  and    inversion  only. 


Ws  =  A(q  +  D  (dh/dt)  -  h  [d  (q+1  )/dt  1  +  vjlllll 

Atq  +  1)  +  (i3h/2) 
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where   h   is   the  inversion  height,   q  and  1  are  vapor  and 

liquid  water  contents,  w'q'  is  the  moisture  flux,  /3  i  s  the 
vertical  gradient  of  (q+i)  in  the  nixed  layer,  and  A-'q-I) 
is  the  difference  between  total  moisture  in  the  nixed  lay^r 


and  total  moisture  immediately  above  thc  mixcd   la; 


In 


reason  *g 


this  evaluation,  the  linear  trend  in  subsidence  during  the 
24  hour  period  pricr  to  an  examination  period  was 
calculated,  and  the  regression  line  applied  at  the  model 
start  time  to  find  the  subsidence  value  used  throughout  the 
predict  ion  peri  od . 

If  the  subsidence  value  was  positive,  then  a  default 
value  of  zero  was  used  to  initialize  the  model.  A  positive 
value  of  subsidence  causes  the  model  to  "blow-up".  The 
his  has  not  been  determined  at  the  tire  of  this 

i 

writ  ing . 

3  •   Sea  Surf  ac_e  Temperature  {S STl 

This  value  was  obtained  from  Acania  data.  The 
SST  at  the  beginning  of  the  period  was  held  constant 
throughout  the  examination  period. 

4 •   ^iHi  Speed 

a.  If  wind  speed  was  not  available  within  3  hours 
of  the  12  hour  point  but  is  available  at  24  hour  point, 
then  a  linear  regression  was  performed  using  the  initial 
and  24  hour  values. 


15 


b.  If  wind  speed  was  not  available  at  the  ?'-  hour 
point,  tut  is  at  the  12  hour  point,  a  linear  repress!:;'1  is 
done  between  the  initial  and  1?  hour  values  and  held 
constant  the  last  12  hours  of  the  period. 

c.  If  wind  speed  was  not  available  within  3  hours 

of  the  1?  or  24  hour  points,  then  the  initial  value  is  held 
constant  throughout. 


ie 


IV.   STATIST  I  CAL/lEIl^  A\ALI3IS  07  THF  MABL  ^:T1 


As  previously  stated,  the  pain  objective  of  this 
thesis   is   to  determine  tre  validity/accuracy  of  trc   YA.EL 

model's   ability   to  predict  ducting  conditions  ever   a  ?A 
hour  period . 

Of  the  4?  cases  used  in  the  verification,  five  cases 
were  omitted  because  during  the  examine t i or  periods  the 
atmosphere  was  stable,  therefore  the  model  physics  do  not 
apply.  Seven  cases  were  omitted  because  tbe  inversior 
height,  Zi,  was  less  than  200  meters,  for  which  ~ase,  the 
model  physics  do  not  apply.    This  left  2?  cases  to  be  used 


t™i 


for  the  model  verification.  These  28  cases  contain  a  total 
of  69  radiosonde  soundings  used  for  verification  purposes. 
Although  the  model  predicts  several  parameters 
(T,q  ,  Zlcl ,  Zi  )  during  the  24  hour  prediction  period,  the 
parameter  under  evaluation  in  this  thesis  is  the 
temperature  inversion  height,  Zi ,  because  this  value  is  the 
height  of  the  top  of  the  atmospheric  duct.  If  this  value 
cannot  be  predicted  accurately,  the  model  is  of  little  use 
to  the  fleet  tactician. 

The   current   method  employed  by   C-eouhysics   Officers 
afloat   to  predict  T  ( tempera  t-ire  )  ,   q  (specific  humidity), 


Zlcl    (lifting   condensation    level), 


an> 


Zi, 


is 


1? 


"Persistence".    Persistence   means  that   nothing  changes. 

For  example,  if  a  radiosonde  launch  at  £"800  indicates  that 
the  inversion  height  is  5  00  meters,  then  a  persistence 
prediction  would  state  that  the  inversion  height  would 
remain  unchanged,  until  the  next  launch  at  which  time  it 
would  he  updated  to  the  current  value. 

K.      STATISTICAL/TEENI  ANALYSIS  DATA  AND  RESULTS 

The  procedures  and  methods  used  for  +hr 
statistical/trend  analysis  of  the  MAEL  model  and 
persistence  will  now  he  presented.  All  of  the  results 
presented  in  this  chapter  will  be  discussed  in  the  chapter 
of  conclus ions . 

1.   Luc t  *rend  Prediction 


$ 


Prior  to  examining  the  inversion  height 
prediction,  it  is  important  to  determine  how  accurately  the 
MA5L  model  and  persistence  predicts  the  trend  of  existing, 
ducts  during  the  24  hour  examination  period.  Predicting 
the  trend  hy  the  model  means,  for  example,  if  an  elevated 
duct  occurs  at  the  initial  sounding  and  then  be^om^s  a 
surface  based  duct,  does  the  model  predict  this  trend? 
Predicting  the  trend  by  persistence  involves  looking  at  thc 
radiosonde  launch  immediately  prior  to  the  initial  sounding 
and  observing  the  trend.  For  example,  if  the  inversion 
height  is  decreasing  during  this  period  then  it  is  assumed, 


IS 


by  persistence,  to  continue  to  decrease  luring  the 
examination  period.  This  segment  of  the  evaluation  does 
not  consider  the  accuracy,  in  meters,  of  the  prediction, 
only  the  trend.  This  means  there  could  he  a  significant 
difference  at  the  verification  time  between  Zi  of  the 
radiosonde  a^d  Zi  of  the  model  or  persistence,  yet  tv-e 
trend  could  have  been  correctly  predicted. 

Considering  all  2.9  cases  under  examination,  the 
model  predicted  the  trend  correctly  60.7%  of  the  tir:e, 
whereas,  persistence  predicted  the  trend  correctly  in  r.~"  .l'\ 
of  the  <~ases  . 

Th*3  28      cases  were  further  subdivided  into   clear 


or   cloudy   sky  conditions,   existing  at  the  on-set 


r  >  p 


examination  period,  to  determine  if  the  initial  sky 
conditions  bad  any  affect  on  the  trend  predictions.  There 
were  18  cloudy  and  iZ  clear  s Ky  cases.  Both  the  model  an: 
persistence  predicted  the  duct  trend  correctly  in  7F.%  (7  of 
1?)  of  the  clear  sky  cases.  Persistence  was  correct  in  5?/: 
(9  of  lc)  of  the  cloudy  sky  cases,  however,  the  model  was 
correct  in  55. ?0'  (12  or  IS)  of  the  cloudy  SKy  cases.  This 
data  is  shown  in  Table  I. 
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TABLE  I 

Duct  Trend  Prediction 

Number.  Type    Model     ^Correct    Pflliiif  nce_ 

28     All     17  correct    c0.7     If  correct 

IS     Cloudy   12  correct    66.7      9  correct 

10     Clear    7  correct    70.0      7  correct 


o,  n   . 


_r r pc t 
57.1 

e  0 . 0 

70.0 


Another  consideration  was,  how  accurately  does 
the  model  and  persistence  predict  sky  conditions  (cloudy, 
^  i  e  a  r  :  at  the  end  of  the  ?  4  hour  period?  T  a  b  1  e  shews 
that  persistence  predicted  the  correct  sky  conditions  in 
50.71  of  the  cases  (17  of  26)  compared  to  45.4%' (12  of  28) 
f  or  the  model . 


TABLE  II 


Prediction  of  Sky  Conditions 
Number    ijO^el     %  Correct    Persistence 
28    13  correct     46.4       17  correct 


%   k2LE£Ct 
60.7 


^  •  l3.IILS.lSS.  Slight  iZil  Prediction 

Table   III  is  a  complete  listing  of  the   raw   data 

used   in   the  analysis.    This  table  gives  the  Late,   Case 

Number,   number   of   hours   since   the   initial  radiosonde 

sounding,    the   radisonde   inversion   height,  "'ALL   model 

prediction  of  Zi,  and  persistence  prediction  of  Zi. 
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TABLE  III 

Complete 

Data  Listing 

for  Zi 

k£l£ 

££2££ 

^o_ur£  Si.nc_g 
Initial 

•   Ei^losonde 

Persis 

25  Sep. 
1975 

2 

12 
25 

1448 
135^ 

1559 
1497 

1466 
14^6 

26  Sep. 
1975 

3 

11 
24 

12°4 
114 

1723 
2065 

1365 
1294 

03  Oct. 
1976 

7 

9 
21 
24 

307 

154 
180 

224 
64 
66 

47^ 
307 

154 

10  Cct. 
1976 

12 

11.5 
24 

193 
336 

193 
142 

250 

193 

11  Oct. 
1975 

11 

11.5 
22.5 

237 

485 

556 
77? 

237 

20  July 
1977 

12 

16.5 
24 

198 
30  5 

747 
671 

510 
196 

21  July 
1Q77 

13 

Is. 5 

23 

24 

190 
546 
339 

48  a 
539 

574 

•  3  0  P, 
190 
54c 

11  May 

1978 

17 

11 

24 

455 
401 
232 

471 

493 

496 

42^ 
455 
401 

14  Pay 
1978 

20 

12.5 
15.5 
24 

176 
1183 
1163 

503 
560 
760 

337 
178 

1163 

19  ^sy 
1978 

21 

12.5 
23.5 

240 
351 

321 
510 

217 
240 

20  May 
1976 

22 

11.5 

24 

465 
628 

545 
51"7 

351 
46  5 

21  May 
1978 

23 

15.5 
24 

685 
566 

695 
1062 

626 
665 
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TABLE    III    (coat.) 


rate 

Case£ 

Houri  Sinc£ 
Initial 

Had  io^onde 

^Ddel 

Persistence 

22  May 
1979" 

24 

4 
10.5 
23.5 

902 
733 

1515 

732 
806 
955 

ess 

902 
733 

01  Aug. 
1979 

25 

14.5 
23.5 

382 
394 

168 
112 

475 

382 

02  Aug. 

1Q78~ 

26 

12 

24.5 

409 
3  50 

424 
461 

384 
409 

03  Aug. 

27 

12 

318 

330 

I07?" 

3  60 

3  4c 

318 

04  Aug. 
1979 

29 

12.5 
24 

299 
241 

50  0 
619 

35? 
299 

05  AU? . 
1979~ 

29 

12 
24.5 

197 
166 

235 
232 

241 
197 

09  Aug. 
1979 

32 

12.5 

19 
24.5 

319 
364 
334 

314 
335 

•  350 

290 
319 

364 

09  Aug. 
1978 

33 

5 
13 

19.5 
24 

339 
343 
445 
154 

288 
254 
234 

220 

334- 

343 
445 

10  Aug. 
1978 

34 

e 

14 

24 

326 
233 
322 

242 
231 
225 

314 
326 
233 

11  Aug. 
1978 

35 

12 
24 

555 
993 

355 
435 

322 
666 

12  Aug. 
1979 

36 

12 
24 

290 
957 

1090 
1436 

045 
28  0 

13  Aug. 
1979 

37 

12 
19 
24.5 

313 
216 
135 

965 
1097 
1233 

867 
313 
216 
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sate 


TASLE  III  (cont. ) 


Cisef  Hours.  Since   EaJJLos_onde  ^odel   Persistence 
Initial 


15  Aug. 
137? 

39 

6 
12 

24 

03  m"ay 
1980 

42 

7 
15 

97   C 

24  ^ay 
1980 

43 

7.5 
13 
23 

05  May 
1980* 

44 

c  m  P 
24.5  . 

792 

254 

30C 

510 

293 

7?? 

767 

345 

510 

482 

532 

544 

292 

73° 

4C2 

506 

?57 

292 

380 

25c 

22  5 

490 

pc.,5 

3  8." 

C  2  "7 

255 

^90 

984 

552 

C,  01  r? 

254 

807 

984 

In  order  to  perform  statistical  analyses  on 
the  predictive  accuracy  cf  the  VA3L  model  ani  persistence, 
the  TTewlet  t -Packard  HP-85  was  used  to  obtair  necessary 
statistical  values.  The  program  "DISTR"  was  used  to 
assemble  input  data  sets  in  ascending  ?rder  a^.r3  to 
calculate  the  mean,  standard  deviation,  and  the  renf?  cf 
the  data  set.  The  program  "Plfs1"  was  used  to  d'ivide  each 
data  set  into  bins  of  equal  width  to  form  a  histogram  of 
the  data  set.  "T^LCT"  was  then  used  to  plot  each  data  set. 
The  program  "PLOTp"  was  used  to  enter  data  pairs,  then 
"TPLOT"  was  used  to  make  a  scatter  plot  of  the  entered 
points.  These  programs  are  listed  in  Appendix  A. 

In  an  effort  to  obtain  more  useful  information, 
five   particular   situations   were   evaluated.     Tfcp   five 
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categories  chosen  were:  1)  all  59  data  points,  cs  57  dats 
points  with  the  high  and  low  values  from  category  1 
excluded,  3)  the  tirre  when  the  verification  was  taken 
during  the  2  4  hour  examination  peri:dt  4)  the  seamen  of  thf 
year  when  the  verificatcn  was  taken,  and  c)  the  initial  sky 
conditions.  Within  categories  3,  4,  and  c  there  were  4, 
3,  and  2  classes,  respectively.  Classes  within  categories 
are  mutually  °rclusive.  Categories  1  and  2  had  only  one 
class  each.  In  all,  there  were  eleven  classes  evaluated, 
each  class  containing  two  data  sets,  one  for  model 


one   for   persistence 


-  j  t 


'he 


ate  s  one 


DISTF.    program. 


associated  classes  are  described  in  Tetle  IV. 

The  Delta  M  and  Delta  F  values  for  each  data  set 
are  listed  in  Appendix:  2,  in  the  output  format  of  the  -?-"" 

Delta  M  and  Delta  F  were  calculated  vy 

subtracting   the   model/persistence  prediction   fro^  vhe 

radiosonde  value  at  that  time.  The  histograms  of  the 
prediction  errors  are  in  Appendix  C. 

Table   V   is   a   compilation   of    the    data  set 

filenames,  the  number  of  data  points,  N,  in  the  set,  and 
the  mean  and  standard  deviation  of  the  set  in  meters. 

Ideally,    the    desired     error    between  the 

radiosonde  sounding  and  the  prediction  value  (model  cr 
persistence)   is  z^ro.    So,   the  desired  mean  of  each  data 

set   is   zero.    The  desired   standard   deviation,   of  the 
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prediction  error,  is,  ideally,  plus  or  minus  *cZ  peters 
(plus  or  minus  120  meters  will  also  be  examined )  .  Any 
inversion  heights  which  cannot  he  predicted  within  plus  or 
minus  5?  meters  is  of  little  or  no  value  to  a  tactician. 
This  predicted  value  is  critical  because  frequencies  which 
are  trapped  within  a  duct  are  a  function  of  duct  thickness 
(Kerr  1951)  and  the  transmitting/receiving  antenna  must  he 
located  within  the  duct  for  maximum  exploitation. 
Therefore,  the  standard  deviation  is  an  important  statistic 
in  that  it  is  a  measure  of  the  precision  of  tee  prediction. 
In  other  words,  it  reveals  hew  much  fluctuation,  a tout  the 
mean,  occurs  for  a  given  data  set.  The  larger  the  standard 
deviation,  the  less  precise  the  prediction. 

The  given  mean  and  standard  deviations  do  not 
clearly  reveal-  the  predictive  worth  of  either  the  model  or 
persistence.  In  order  to  make  more  concise  statements 
ahout  the  two  predictive  methods,  it  will  he  useful  to 
perform  hypotheses  tests  utilizing  this  data  and  a  level  of 
significance ,  a   ,  of  0.01. 

a.   Hypothesis  Concerning  Cne  fean 

Hypothesis  testing  is  done  to  determine  whether  or 
not  a  given  statement  concerning  a  data  set  can  re  accepted 
with  a  probability  a  of  risking  a  Type  I  error  [Ref .  3:  p. 
195].   This  is  referred  to  as  the  null  hypothesis,  Ho. 


CO 


TAILE    IV 


Ca tegory /Class    Description 


Category   ^uniher 
1 


Glass    Nurrfcer 

1 
2 


13 
11 


pescripti on 

All  59  data  points  usee 

67  data  point  s  uscd  . 
The  high  and  lev 
values  from  data  set  1 
were  <=xciudei. 

Radiosonde  soundings 
whi  ch  occurred  4  to 
10.99  hours  after 
initial  launch 

Radiosonde  soundings 
which  occurred  11  to 

13  hours  after 
initial  launch 

Radiosonde  sounding? 
which  occurred  l77 .71 
to  22.39  hours  after 
initial  launch 

Radiosonde  soundings 
which  occurred  23  to  2c 
hours  after  initial 
launch 

Radiosonde  soundings 
taken  during  March, 
April,  and  May 

Radiosonde  soundings 
taken  du ring  June , 
July,  and  August 

Radiosonde  soundings 
taken  during  September, 
October,  and  Never her 

Sky  conditions  clear 

Sky  conditions  cloudy 


TAELE  '/ 


yean  end  Standard  Deviation 

Standard  D  e  v  i  a  t  i 


£ita  Set 

N 

Mean 

Filename 

vOQ^ll 

59 

S0.S12 

?ers  is  tl 

69 

-4  .555 

yodel2 

67 

62.806 

?ersi  st2 

67 

-11  .581 

vodel3 

10 

-121  .100 

Persist3 

10 

-57.500 

!*odel4 

19 

149.056 

Persists 

IS 

93.275 

Models 

12 

158.917 

Persists 

12 

40.917 

Model? 

27 

30  .704 

Persist 6 

27 

-58.519 

Pod»17 

22 

55  .  554 

Pers  ist7 

22 

-46.591 

Mode 18 

35 

51.085 

Persists 

35 

0.000 

Model9 

10 

81.900 

Persist9 

10 

27.300 

vodell0 

25 

22.920 

?ersistl0 

25 

-38.720 

Modelll 

44 

140.023 

Persistll 

44 

7.670 

410  .02* 


425 
242 

091 

947 


335 

217 

■i    /~\  /-* 
XJC 

183 

261 

209 


145. 5^8 


353 
223, 

1  CC 

.  9 '.-?  5 

287 
203 

.037 

390, 

241 

.  0  52 
.9  52 

198. 
124, 

.228 

.982 

450 
335, 

.58  0 
.523 

372 
248, 

777 

.5  00 
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Figure  4.   Type  I  Error 

Given  the  distribution  in  Figure  4,  a  Type  I 
error  would  be  the  probability  a  of  rejecting  the 
hypothesis  wren  it  is  true. 

The  null  hypothesis  chosen  was,  Ho:  \i  -  f*o, 
which  in  words  states,  is  the  mean  of  Delta  !VTelta  ?  equal 
to  mo,  a  specified  value  for  the  desired  mean.  Since 
Delta  M  and  Delta  P  are  the  errors  associated  with  the 
differcnce  between  the  model  or  persistence  prediction  ~nd 
the  radiosonde  value  and  the  desired  error  to  be  zero,  then 
Mo  was  chosen  to  equal  zero.  Another  way  of  stating 
this-  would  be:  Can  it  be  said  with  reasonable  confidence 
(99^)  that  the  mean  of  the  error  for  each  predictive  method 
be  assumed  to  he  zero? 

In  order  to  carry  out  this  test,  the  Student's 
t-test  was  used.  The  assumptions  required  for  this  test 
are  that  the  actual  standard  deviation,  a  ,  be  unknown  and 
that  the  sample  comes  from  a  normal  population.  The  t 
stat  ist  ic  is  , 
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t  = 


X  -   MO 


SZ/n 


where  x  equals  the  mean  of  the  sample,  \i  o=  0.0,  S  is  the 
sample  standard  deviation,  and  n  is  the  number  of  data 
points  in  the  set.  This  calculated  value  of  t  was  compared 
to  a  value  obtained  from  a  t-statistic  table.  This  value 
was  obtained  by  entering  the  table  with  the  number  of 
degrees  of  freedom,  in  this  case  n-1,  and  a  /2  =  0.005 
because  this  is  a  two-tail  test.  If  the  absolute  value  of 
t -calculated  is  greater  than  t -table  than  the  hypothesis  is 
accepted.  Figure  5  shows  the  acceptance/rejection  regions 
for  a  two-tail  test  [Pef.  3:   pp.  211-214]. 

The  results  of  the  hypothesis  test,  Ho:  m  =  mo   for 
each  data  set  are  shown  in  Table  VI. 


Reject 

Accept 

Reject 

• 

/ '  y///'  ■ 

/  '-W.'/Zy 

- . .  .  ■  .  , .. 

v. 

/       \ 

/ ' '  ■  ■  \ 

\ 

V//A/,- 

-   

Figure  5.   Acceptance/Rejection  Regions 
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TAB 

LI  VI 

Fo:  M  = 

fJ-0     =     0.0 

Filename    n 

trcalc 

t-ta>ie 

Modell 

69 

1.540 

2.576 

Persistl 

59 

-0.135 

-2.575 

Model? 

57 

1.533 

2.576 

Persist2 

57 

-0.44S 

-2.576 

Mod<=13 

1? 

-1  .930 

-3.2  50 

Persists 

10 

-0.990 

-3.250 

Model4 

ie 

2.415 

2  .  898 

Persist-1 

18 

i.ese 

2.898 

vodel5 

12 

1.552 

3.105 

Pers  is t5 

12 

0.952 

3.106 

Model6 

27 

0.415 

2.779 

Pers  is  to 

27 

-2.053 

-2.779 

Model? 

22 

0.926 

2.831 

Persist? 

22 

-1.123 

-2.831 

vodel8 

35 

0.927 

2.576 

Persists 

35 

0.000 

2.576 

Model9 

10 

1.307 

3.250 

Persist  9 

10 

0.691 

3 . 2  50 

ModellS 

25 

0.254 

2.797 

Persistl0 

25 

-0  .577 

-2.797 

Modelll 

44 

2.490 

2.576 

Persistll 

-14 

0.210 

2.576 

Acc£ttZPp jeci 

A ccept 
Accept 

Accept 
Accept 

Accept 

£  p  n  c   T-.  t 

A  c.c  e  p  t 
Accent 


Accept 
Accent 


Accept 
Accept 

Accept 
Accept 

Accept 
Accept 

Accept 
Accept 

Accept 
Accept 

Accept 
Accept 
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t.   Hypothesis  Concerning  One  Variance 
The  null  hypothesis   considered   in   this   section 
will   be  tc  determine  if  the  population  standard  deviation 
is  equal  tc  a  specified  constant  (Ho:  o  -     or,).        This  test 
is  conducted  by  calculating  the  Cni-Square  statistic, 


i      (n-l)S 

*    =  2 

oo 

where   n   is   the   numher  of   data   ooints   in   the   samtle 

population,     is  the  variance  of  the  sample,   and  oq      is 

the  specified  ''desired)  variance  (variance  is  the  square  of 

the   standard  deviation).    The  Chi-Square  calculated  value 

is   then   compared  to  a  Chi-Square  value   obtained   from   a 

table  in  Ref  .  3,  using  n-1  degrees  of  freedom  and  a=  0.01. 

If  the  calculated  value  is  less  than  X afe       0T   fcreat°r  than 

v  ,     ,   the   null  hypothesis  is  accepted.    The   desired 
i-o/l 

standard  deviation  was  chosen  to  be  o  o  =•  52  meters, 
however,  ao  ~  l%2  meters  was  also  tested.  The  results  of 
these  tests  are  given  in  Table  VII  [Eef.  3:   pp.  233-235]. 

c.   Hypothesis  Concerning  Two  Means 

In  order  to  determine  if  there  is  a  significant 
difference  between  the  two  means,  the  model  and 
persistence,  the  null  hypothesis,  Ho:  ^  =  ftp  must  he 
examined.    Along  with   the   assumption   of   normality,   a 

statement   must  be  made  about  the  variances  cf  the  two  data 

Z 
sets.    Namely,   can   it  be  assumed  that  the  variances   am 
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TABLT  VII 


H  o  :  o    =  a  o 


ao  =  5?  and  10?  meters 


Rilena^e. 

Mode 11 

-ers  i  stl 
vcdeI2 
Persist 2 
Mode  13 
Persists 
Model4 
Persists 
Modelc 

Persi  s:5 

vodel5 

Persist? 

vodel7 

Persist? 

Models 

Persists 

Model9 

?ersist9 

M  od  e  1 1 0 


1.7x10 
i 

7.8x10 

1.1x10 

4.7x10 

3.9x10 

3.4x10 

6.9x10 
i 

4.4x10 
t 

1.1x10 
2.2x10 

c 

1.5x10 

i 

5.0x10 
i 

8.3x10 
i 

4. 1x10 

c 

1.5x10 
t 

5.9x10 
/ 

3.9x10 
1.6x10 

c 

2.0x10 


5 


z 
X-cal^ 

50 

ulated 
102 

liable 

4570 

1140 

102.3 

2130 

532 

102.3 

2980 

746 

99.5 

1250 

312 

99.  5 

141 

7C, 

p7  a 

122 

30 

23.5 

466 

H7 

35.7 

299 

74 

TR  7 

4S8 

122 

26. 8 

97 

24 

26  .8 

1530 

332 

48.3 

522 

131 

48.3 

693 

173 

41.4 

346 

S* 

41.4 

2070 

517 

56.1 

796 

199 

56.  1 

141 

35 

23.6 

55 

14 

23.6 

1950 

487 

45.6 

£££££i/E£JJ£rt 

50  100 

Reject  Reject 

Reject  Reject 

Reject  ^  e  j  e  c  t 

Feject  Reject 

Reject  Re jec t 

Reject  P.  e  i  e  c  t 

Reject  Reject 

Reject  Reject 

Reject  Reject 

Feject  Accept 

Reject  Reject 

Reject  Reject 

Reject  Reject 

Reject  Reject 

Reject  Reject 

Reject  Reject 

Feject  Reject 

Reject  Accept 

Reject  Reject 
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TOLE    VII    (cont  .  ) 


Persists    l.lxl? 


tf ode  111        1.4x1? 


Persist  11  6.2x10 


X-calcu] 
50 

.ated 
100 

X-iable 

ies0 

271 

45.5 

2390 

599 

71.0 

1060 

2c5 

71  .0 

50   100" 
Reject   Reject 
Reject   Reject 
Reject   Reject 


i  i        2 

and   ^p   are  equal?    The  null  hypothesis  Ho:    CT  m  =   ffp 

must  he  completed  to  answer  this.  If  it  can  be  said  that 
they  are  equal,  (i.e.  accept  the  hypothesis),  then  a  two- 
sample  t  test  can  be  used  to  test  the  means.  If  the 
hypothesis  that  the  variances  are  equal  is  rejected,  then  a 
paired-sample  t  test  must  be  used  to  test  the  means. 

In  order  to  perform  the   hypothesis   test,   an   F- 
statistic  was  calculated  by, 


F  = 


1 

Sm 


where  Sm  is  the  variance  of  the  model  data  and  Sp  is  the 
variance  of  the  persistence  data.  Using  a=  0.01,  if  F- 
calc  was  less  than  F-table  then  the  hypothesis.  Ho:  am  = 
ap  must  be  accepted.  The  results  are  shown  in  Table 
VIII. 
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TABLE  VIII 


He: 

1 

O  IT 

2 
=       *P 

Vata    Set    £ 

IlSalCUlated 

F-tatle 

1 

2.1? 

1.66 

2 

2 .  39 

1.66 

3 

1.16 

5.3  5 

4 

.1.56 

2  .c2 

c 

5.00 

4.54 

a 

2.92 

2.52 

7 

2.00 

2.88 

c 

2.50 

2.30 

a 

2.52 

5.3^ 

1? 

1.80 

2.65 

11 

2.25 

2.0  5 

Accept /Eejfj 
r'  £  j  e  c  t 
F  e  ,1  e  c  t 
Accept 
Accept 
Reject 
F,e  jf  r  t 
Accept 
F  e  j  e  c  t 
Accept 

Reject 


7  a 


five      data      sets      a  c  c  e  j. 


variances      were      then      tested      with      the      hypothesis,         r~o: 
H  m     -     /ip      =0.0.        The      two-sample    test    was   use:      where, 
after    reduction    of    the   general    formula     [Eef.    3:       p.    21?], 


(  M  p-   -     MP    )      n  (  n -1  ) 


t    = 


(n-1)  [Srn 


ft   i 


with   a  =  0.01  snd   the  decrees  of  freedom  being 
The  results  are  given  in  Table  IX. 
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TABLE  IX 


Ho:  V-t    -     /x p  =  0 . 2 


£ata  Set* 
3 
4 

7 
o 

10 


A             A 

Si   ±    St 

lioalc. 

tztab.1 

4. 

7  .  3x  1 0 

-0.71 

_o     2  O 

c. 

1.1x10 

0.59 

2.5- 

c 

1.2x10 

1.37 

2.56 

4 

5.5x10 

0.70 

2.92 

5 

3.2x10 

0.22 

2.58 

Acce^l/Ili  jec_t 
Accept 
.a,  r  c  e  p  t 
Accept 
Accept 
Accept 


The  six  data  sets  rejected  ?or  not  being  able  to 
make  the  assumption  that  the  variances  are  equal  had  their 
means  tested  for  equality  using  the  paired-sample  t  test. 
In  order  to  conduct  this  test,  a  data  set  containing  the 
values,   Delta   M   -  Telta  P,   was  compiled  and  a  mean   a nd 
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standard  deviation  for  each  data  sct  was  calculated  using 
"riSTR".  The  desired  rear,  |io,  cf  each  date  set  is  once 
again  zero.  The  t  statistic  was  calculated  for  e:c'r.  of  the 
6  sets  as  in  Section  A.  The  six  paired  data  sets  are 
listed  in  Appendix  r  ard  the  results  of  the  hypothesis, 
Ho:  \i   =    no   =  0.0,  are  listed  in  Table  X. 

TAPLP  X 


£ata  s°i  £ 
1 

2 

5 

c 

11 


F  o  :  M  ■ 

izc^lcula t e; 
2.91 
3.31 
1 .55 
2 .  59 
1.54 


Mo  =  0.0 
trtatlf 

2.55 
2  . cc 
3.06 
2.75 
2  .  5£ 
2.56 


A  r  c  e 7  ■  t  /  ?.  e  ;  * 


Reject 
Reject 
Accept 
A.  c  c  e  p  t 
A  c  c  e  p  t 
Reject 


d.  Prediction  Error  vs.  Verification  Time 
Figures  5  and  7  are  plots  of  the  Prediction  Error 
vs.  Verification  Time.  With  zero  "being  the  desired  error, 
it  is  clear  that  as  the  model  progresses  further  into  the 
24  hour  prediction  period,  that  the  prediction  error 
increases  significantly.  The  standard  deviation  varies 
from   199.091    meters    to   363.166    meters    (from    Table   V). 
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The  Persistence  Prediction  Error  plot  reveals 
much  less  scattering  of  the  data  points,  while  there  is  no 
obvious  trend  to  the  prediction  error. 

e.  Prediction  Error  vs.  Month  of  Observation 
Figures  8  and  9  are  plots  of  the  69  data  points 
relative  to  the  month  in  which  they  ere  taken.  The  /-axis 
of  these  two  figures  are  labeled  by  numbers,  where  1 
corresponds  to  January  through  12  which  is  December.  It  ii 
obvious  that  there  is  much  mere  scatter  associated  with  tve 
model's  prediction  error  during  the  summer  months  than  at 
other  times. 
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V.   CONCLUSIONS 

A.  DUCT  TREND  PREDICTION 

From  a  statistical  viewpoint  there  is  no  significant 
difference  between  the  MABL  model  and  persistence  in 
predicting  the  duct  trend  or  sky  (clouds/no  clours) 
conditions.  There  are  seme  differences  in  their  correct 
prediction  percentages,  however,  net  enough  to  cheese  one 
method  ever  the  other.  The  results  clearly  indicate  that, 
given  sky  conditions  are  clear  at  the  initial  start  time, 
the  duct  trend  can  be  accurately  predicted  by  either 
method.  Cloud  formation  appears  to  cause  problems  in  the 
model's  predictive  ability  and  should  be  further 
invest  igat  ed  . 

B.  HYPOTHESES  TESTS 

Although  the  hypothesis  test  of  Section  Ilia.  shewed 
that  the  mean  of  each  data  set  had  to  be  acceptea  as  being 
equal  to  zero,  this  can  prove  to  be  misleading.  This 
result  should  be  easily  expected  considering  the  large 
standard  deviations  in  the  data.  The  hypothesis  test  which 
contains  extremely  important  results,  from  the  tacticians 
viewpoint,  is  the  test  of  Section  Illb.  Remembering,  that 
standard   deviation   is  a  measure  of  the  precision  .of   the 


42 


prediction,  it  car.  he  stated  quite  assuredly  that  the 
inversion  height  cannot  "be  predicted  accurately,  eithcr  by 
the  model  or  persistence,  to  the  degree  required  for 
tactical  employment.  In  all  cases,  the  accuracy  of 
persistence  was  tetter  than  the  model,  hut  still  not  good 
enough  for  tactical  use. 

The  hypothesis  tests  concerning  the  equality  of  the 
means  revealed  that  there  is  no  significant  difference 
between  the  populations.  Once  again,  this  can  he 
misleading  due  to  the  large  standard  deviations 
enc  oun tered . 


C.   PREDICTION  ERROR  VS.  VERIFICATION  TIME 

Figure  5  clearly  indicates  the  model's  predictive 
atility  decreases  significantly  the  further  it  progresses 
into  the  24  hour  period.  Figure  7  shows  much  less  scatter 
associated  with  the  persistence  predictions.  These  plots 
seem  to  indicate  that  the  model  was  expecting/predicting  a 
much  greater  change  in  Zi  than  actually  occurred. 
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T.   PREDICTION  FRROR  VS.  MONTH  OF  OBSERVATION 

Figure  6  indicates  the  model's  predictive  ability  is 
significantly  reduced  during  the  summer  months.  There  is 
much  more  scatter  associated  with  the  summer  months  than 
others.  The  month  of  observation  had  no  discernible  affect 
on  persistence  as  seen  in  Figure  9. 
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E.   SKY  CONDITIONS 

Initial  sky  conditions  had  no  apparent  affect  en  the 
nodel  or  persistence. 

?.   GENERAL  CONCLUSIONS 

The  results  appear  to  point  toward  a  conclusion  that 
the  MA  EL  model  and  persistence  are  not  useable,  from  a 
tactical  viewpoint.  More  testing  and  statistical  work 
should  be  done  before  a  judgement  can  be  made  about  the 
model  physics.  The  next  chapter  deals  with  approaches 
which  should  be  taken  in  an  effort  to  make  a  more  valid 
evaluation  of  the  MAP-L  model. 


t 
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VI.   RECOMMENDATIONS 

Proper  experimental  planning  can  give  a  reasonable 
assurance  that  the  results  of  an  experiment  will  provide 
clear-cut  answers  to  questions  under  investigation. 

The  examination  in  this  thesis  made  the  major 
assumption  that  the  radiosonde  data  was  "perfect"  and  that 
any  results  would  indicate  the  precision  of  the  models' 
predictive  ability.  However,  it  is  widely  known,  and 
accepted,  that  radiosonde  data  are  not  "perfect"  because 
the  radiosone  is  not  sensitive  enough  to  abrupt  changes  in 
atmospheric  conditions  as  it  ascends.  In  order  tc 
determine  if  the  error  contributed  by  the  radiosonde  is 
significant  and  if  the  radiosonde  data  differs 
significantly  from  an  alternate  method  of  collecting  the 
same  type  of  data,  an  Analysis  of  Variance  (ANOVA)  study 
should  be  done  involving  three  separate  comparisons.  A 
complete  discussion  of  this  method  can  be  found  in  Chapter 
12  of  Fef.  3.  The  mathematics  of  this  method  will  not  le 
discussed  here,  but  suggested  procedures  for  conducting  the 
necessary  data  collection  will  be  elaborated  upon. 

The  three  ANOVA  comparisons  recommended  are  delineated 
in  the  following  subsections. 
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A.  radiosonde  vs.  radiosonde 

To  determine  if  the  data  received  by  radiosondes  is 
consistent,  a  test  should  be  conducted  in  which  two 
radiosondes  are  launched  within  minutes  of  each  other  at 
each  observation  time.  Ensuring  that  both  radiosondes  are 
identically  calibrated  and  launched  from  a  weather  deck 
where  there  will  be  little  possibility  of  collision  with 
other  structures  during  ascent,  the  first  radiosonde  should 
be  launched  and  then  the  second  launched  a  few  minutes 
later.  Eual  launches  should  be  made  at  each  measurement 
period.  The  number  of  dual  launches  made  should  be  limited 
to  between  5  and  12  because  of  cost  and  a  large  number  of 
dual  launches  will  not  be  necessary  to  check  f 
consistency  in  data.  It  is  imperative  that  the  ANOVA  tests 
be  carried  out  even  if  by  observa  ion  the  data  appears 
consistent.    This   data  can  then  be  used  to   evaluate   the 
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consistency   in   received    radiosonde   data. 


Regardless  whether  the  tests  determine  if  the  data  from  two 
different  radiosondes  is  consistent  or  net  consistent, 
further  tests  will  be  necessary. 

E.   RFERACTOMETER  VS.  REERACTOMETER 

The  AN/AMF-3  (XAN-5)  Airborne  Microwave  Refract ometer 
(A!*R)  is  currently  used  in  some  Navy  aircraft  to  ottain  the 
same   type   of   information   which   radiosondes   provide. 
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Simultaneous  data  should  be  collected  by  two  aircraft  tc 
determine  the  degree  of  consistency  in  refract  ore ter 
measurements .  This  test  will  require  the  same  aircraft  tc 
carry  two  different  refrectometers  or  two  aircraft  each 
carrying  an  AVR  and  launched  within  minutes  of  each  other 
and  flying  the  same  pattern.  As  a  result  of  the  ANCVA  test 
of  this  data,  if  the  consistency  is  good  then  the 
Ref Tactometer  can  he  used  in  a  test  comparing  its'  data  to 
the  radiosonde. 

C.  RADIOSONDE  VS.  REFRACTOMETIR 

In  order  to  determine  if  there  is  any  significant 
differcnce  between  radiosonde  data  and  ref Tactometer  data, 
data  should  he  collected  A thin  32-6?  minutes  of  each 
other.  If  there  is  significant  difference  in  the  data  then 
careful  consideration  should  he  taken  to  determine  which, 
if  either,  method  provides  the  more  accurate  data. 

Once  these  tests  have  been  done  and  the  statistical 
results  determined,  then  much  clearer  statements  can  be 
made  about  the  model's  predictive  ability  because  more  will 
be  known  about  the  quality  of  the  input  data. 

D.  GENERAL  RECOMMENDATIONS 

The  other  parameters  which  are  predicted  by  the  model 
should  undergo  similar  statistical  analysis  as  that  which 
has   been   done   on  the  parameter   Zi.    This   analysis   is 
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necessary  so  that  a  clearer  understanding  of  the  model's 
predictive  abilities  is  icnown. 

Some  caution  and  consideration  should  be  given  to  the 
results  of  the  three  tests  suggested  above,  because  the 
results  for  any  or  all  of  the  tests  for  one  geographical 
location  and/or  tine  of  year  of  the  test  may  agree  or 
disagree  dees  not  necessarily  mean  the  results  will  or  will 
not  agree  for  another  location  or  different  tire. 

No  attempt  should  he  made  to  alter  the  VAEL  model 
until  the  ANOVA  tests  are  done.  A  qualitative  and 
quantitative  assessment  of  the  input  data  is  imperative 
before  a  fair  evaluation  of  the  MABL  model  can  be  nade. 
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APPENDIX  A 


FP-S5  PROGRAMS 


The  programs  used  in  this  thesis  are  listed  on  the 
following  pages.  Each  program  can  be  used  on  an  FP-E5  by 
first  starting  the  machine  and  placing  a  disk  with  the 
programs  on  it  into  the  disk  drive.  Next  type,  """ASS 
STORAGE  IS  ":D?03",  then  hit  the  END  LINE  key,  then  type, 
LOAD  "program  name"  and  then  the  END  LINE  key.  The  END  LINE 
key  muse  be  used  after  an  entry.  Once  the  program  is 
loaded,  hit  the  RUN  key  and  the  program  will  prompt  for 
inputs.  ^ 

The  program  DISTR  will  first  prompt  for  a  data 
filename.  Upon  entering  the  filename,  the  next  prompt  will 
be  to  enter  the  number  of  data  points.  The  data  can  be 
entered  in  any  order.  Once  the  last  point  has  been  entered 
the  program  ranks  the  data  in  ascending  order,  computes  the 
mean,  standard  deviation,  maximum  and  minimum  values  of  the 
data,  and  the  range.  It  then  prints  the  data  plus  the 
results  of  the  computations. 

The  program  "BIN"  is  used  to  divide  the  data  into  bins 
of  equal  width  in  order  to  make  a  histogram  plot.  The 
program  prompts  for  the  data  filename,  then  the 
distribution  filename.   The  data  and  distribution  filenames 
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must  differ.  The  next  pronpt  is  for  the  number  of  tins. 
This  value  is  any  integer  value  less  than  the  number  of 
data  points.  The  number  of  "bins  chosen  should  he  such  that 
there  are  at  least  two  data  points  per  bin  whenever 
possible . 

Upon  completion  of  the  "BIN"  program,  "TPLOT"  can  be 
loaded  for  use  in  making  a  histogram  plot.  The  first- 
prompt  is  for  the  data  filename,  however,  enter  the 
distribution  filename  for  which  a  plot  is  desired.  P 
prompt  for  the  title  and  subtitle  of  the  plot  are  displayed 
next.  If  either  or  both  of  these  are  not  required,  hit  the 
ENE  LINE  key  and  the  next  prompt  will  be  displayed.  The 
next  entries  are  XM  LXJ^AX  .  YM I N ,  YMAX .  These  are  the 
maximum  and  minumum  valuec  for  the  horizontal  and  vertical 
axis,  respectively.  The  four  values  should  be  entered  on 
the  same  line  separated  by  commas  before  keying  END  LINE 
The  next  prompt  is  X  AXIS  LOG  SCALE  ENTEE  Y  or  N.  Upon 
answering  yes  or  no,  the  next  prompt  is 
XINTV,YINTV,XINTC,YINTC .  These  values  are  the  x  and  y 
intervals  for  labeling  of  the  x  and  y  axis  and  the  x  and  y 
axis  intercepts.  The  *  and  y  axis  titles  are  requested 
next.  Upon  completing  these  entries,  a  plot  is  made  on  the 
screen.  When  the  plot  is  done  hit  the  CONTINUE  key  and  the 
number  p  'for  a  hard  copy)  and  it  will  plot  on  the  plotter. 
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The  PLOT?  program  should,  be  use-*  if  x  and  y  paired 
data  is  required  to  be  plotted.  The  first  prompt  is  for 
the  number  of  paired  data  sets.  The  next  q^iery  is  for  the 
data  filename.  The  next  pronpt  will  he  for  the  paired  data 
sets.    Fnter   the  x  coordinate  first  followed  by   a   comma 


then   the  v  coordinate  followed  by  the  EM  LINE  fcev . 


7Ti  rv 


completion   of  entering  the  data  sets,   "T?L0T"  can  he  usee 
to  make  a  scatter  plot. 
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B  I  N 


10  DIM  DC250) .FC250) 

20  DISP  'DATA  FILENAME" 

30  INPUT  D* 

40  DISP  "DISTRIBUTION  FILENAME" 

50  INPUT  F% 

60  PRINT  "DISTRIBUTION  FILENAME 

'  ,F$ 

70  DISP  "NO  OF  BINS" 

80  INPUT  G 

90  ASSIGN/  1  TO  Dt 

100  READ/  1  i  DCO) 

1 10  FOR  1-1  TO  DCO) 

120  READ/  1  ,  0(!) 

130  NEXT  I 

140  B«CDCO)«l)  DIV  8*1 

150  8-8*8 

160  CREATE  FI.16.B 
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210  PRINT  "BIN  PT,  NO.,  LIMIT,  D 

.  PT.  " 
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350  END 
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140  IF  P-4  THEN  200  !  MORE  DATA 

SAME  PLOT 

150  IF  P-5  THEN 

160  GOSUB  380  » 

170  IF  P-2  THEN 


200  !  MULTIPLOT 

TITLES 

190  !  CHG  TITLES 


180 
190 
200 
210 
220 
230 
240 
250 
260 
270 

280 
290 
300 
310 
320 
330 
340 


X/Y  AXES 
DRAW  AXES 
•  PLOT  DATA 


GOSUB  440 
GOSUB  770 
GOSUB  1430 
PAUSE 

1  SELECT  PLOT 
OISP  'CHOOSE  PARAMETER' 
1-REDRAW' 


DISP  ' 
DISP  ' 
DISP  ' 
DISP  " 

DISP  • 
DISP  ' 
DISP  ' 
DISP  ' 
INPUT  P 
IF  P-4  THEN 
IF  P-8  THEN 


2-CHG  TITLES' 
3-CHG  X/Y  AXES' 
4-MORE  DATA  SAME 

PLOT' 
5-MULTIPLOT' 
6-NEW  PLOT' 
7-END  PLOTTING* 
8-HARD  COPY' 


90 
PLOTTER 


350  IF  P-l  THEN  PLOTTER  IS  705  • 

GOTO  80 
360  IF  P-7  THEN  1680  !  END 

GOTO  70 

!  TITLES 

DISP  'TITLE' 

INPUT  A* 

DISP  'SUBTITLE' 

INPUT  Alt 

RETURN 

i  X/Y  AXE S 

DISP  'INPUT  XM1N,XMAX,YMIN,Y 


370 
380 
390 
400 
410 
420 
430 
440 
450 
MAX 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 


INPUT  XI  ,X2,Y1 ,Y2 

DISP  'X  AXIS  LOG  SCALE?' 
DISP  'ENTER  Y  OR  N' 
INPUT  QOt 

THEN  560 


IS  705  • 


IF  QOiO'Y' 

X3-0 

X4-X1 

DISP  'YINTV.Y1NTC 

INPUT  Y3.Y4 

GOTO  660 

DISP  'INPUT  XINTV.YINTV.XINT 
C.YINTC" 

570  INPUT  X3.Y3.X4.Y4 
580  IF  P-0  THEN  620 
590  DISP  'CHG  AXIS  LABELS?  Y/N' 
600  INPUT  01$ 
610  IF  Qlt«'N'  THEN  660 
620  DISP  'XAXIS  TITLE" 
630  INPUT  Bt 
640  DISP  'YAXIS  TITLE' 
650  INPUT  Ct 
660  RETURN 
670  '  READ  DATA 
680  DISP  'DATA  FILENAME' 
690  INPUT  Ft 
700  ASSIGN*  1  TO  Ft 
710  READ/  1  j  N9 
720  N9-2«N9 
730  FOR  1-1  TO  N9 


GOTO  80 
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T  P  L  0  T 


(coat.) 


740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

26, 

860 

870 

880 

890 

900 

-2 

910 

920 

930 

940 

95  0 

960 

970 

-2 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 


READ*  1  |  DCI) 

NEXT  I 

RETURN 

!  DRAW  AXES 

LOCATE  30,110,32,89 

C5IZE  3 

'  DRAW  LINEAR'' 

IF  QOtO'Y"  THEN  GOSUB  1010 

»  OR  DRAW  SEMILOG 

IF  QOt-'Y*  THEN  GOSUB  1160 

!  PUT  LABELS 

PLOT  <X1*X2)/2,Y2»3*CY2-Y1)/ 

2 

LORG  6 

CSIZE  4 

LABEL  A* 

LABEL  Ai« 

PLOT  (Xl«X2)/2,Yl-(Y2-Yl)/8, 


LORG  4 

LABEL  B* 

DEG 

LDIR 

LORG 

LORG 

PLOT 


90 

6 

5 

X1-CX2 


Xl)/8,(Yl*Y2)/2, 


LABEL  C« 
LDIR  0 

RETURN 

»  DRAW  LINEAR  AXIES 

SCALE  X1,X2,Y1,Y2 

AXES  X3,Y3,X4,Y4 

LORG  5 

FOR  X-Xl  TO  X2  STEP  X3 

PLOT  X.Y1-CY2-YD/55 


LABEL  USING 
NEXT  X 
LORG  8 
FOR  Y-Yl  TO 
PLOT  X1-CX2- 
LABEL  USING 
NEXT  Y 


J 


Y2    STEP    Y3 
XD/40.Y 
'K'     |     Y 


1140 

RETURN 

1150 

>  DRAW  L0GX-AX1S 

1160 

LORG  6 

1170 

IF  Lt-'SET*  THEN  1210 

1180 

Lt-'SET" 

1190 

X1-LGTCX1) 

1200 

X2-LGTCX2) 

1210 

FOR  1-0  TO  8 

1220 

K1CI)-LGTCI*2) 

1230 

NEXT  I 

1240 

SCALE  X1.X2.Y1.Y2 

1250 

MOVE  XI, Yl 

1260 

DRAW  X2.Y1 

1270 

FOR  I-Xl  TO  X2 

1280 

MOVE  I  .Y1-CY2-YD/55 

1290 

LABEL  USING  mK'     ,  10*1 

1300 

FOR  J-0  TO  8 

131  0 

K2-K1(J)»I 

1320 

MOVE  K2.Y1 

1330 

DRAW  K2,Y1*(Y2-Y1)/100 

1340 

NEXT  J 

135  0 

NEXT  I 

1360 

YAXIS  XI, Y3 

1370 

LORG  8 

1380 

FOR  Y-Yl  TO  Y2  STEP  Y3 

1390 

PLOT  X1-CX2-X1)/40,Y 

1400 

LABEL  USING  'K'  {  Y 

1410 

NEXT  Y 

1420 

RETURN 

1430 

!  PLOT  DATA 

1440 

DISP  'L1NETYPE  ?" 

1450 

DISP  •   1  SOLID' 

1460 

DISP  '   2  END  PT  ONLY' 

1470 

DISP  '   3  DOTS' 

1480 

DISP  '   4  SDASH' 

1490 

DISP  '   5  LDASH' 

1500 

DISP  '   6  DASH  DOT' 

1510 

DISP  '   7  LDASH  SDASH' 

1520 

DISP  '   8  LDASH  SDASH  SDASH 

1530    INPUT   L 


55 


T  P  L  0  T    (cont.  ) 


15*0  L1NETYPE  L 

1550  IF  Q0*O'Y'  THEN  1580 

1560  MOVE  LGTCD(D)  ,D(2) 

1570  GOTO  1590 

1580  MOVE  DC1) ,DC2) 

1590  FOR  1-1  TO  H9  STEP  2 

1600"  IF  QOtO'Y"  THEH  1630 

1610  TLOT  LGTCCCm  ,DC1«1> 

1620  GOTO  1640 

1630  PLOT  DCI)  ,DCI«n 

1640  NLXT  I 

1650  PEMUP 

1660  PENUP 

1670  WZJURH 

1680  STOP  ' 

1690  END 
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P  L  0  T  P 


10  DISP  'NO.  OF  POINTS' 

20  INPUT  N 

30  DIM  DC200) 

40  DISP  'DATA  FILENAME' 

50  INPUT  D« 

60  B-N  DIV  8*1 

70  B-B«8 

80  CREATE  Dt,16,B 

90  D(0)-N 

100  FOR  1-1  TO  2»N  STEP 

110  DISP  'INPUT  X,Y' 

120  INPUT  OCI)  ,D(  !♦!) 

130  NEXT  I 

140  ASSIGN/  1  TO  D* 

ISO  FOR  1-0  TO  2*N 

160  PRINT*  1  ,  DCI) 

170  NEXT  I 

180  ASSIGN/  1  TO  I 

190  END 
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APPENDIX  B 
DILTA  M/DELTA  ?  VALUES 

The  Delta  v  and  Delta  ?  values  for  each  data  set  are 
listed  on  the  following  pages  in  the  output  format  of  the 
IIP-35  "DISTr"  program.  The  Delta  M  values  were  attained  "by 
subtracting  the  radiosonde  value  for  Zi  from  the  value 
predicted  by  the  model  at  that  time.  The  Delta  ?  values 
were  obtained  by  subtracting  the  current  radiosonde  value 
from  the  previous  value. 
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DATA    FILENAME   M0DEL1 


603 

-579 

-56  0 

548 

-538 

-422 

401 

-272 

252 

-234 

-234 

214 

-211 

-211 

194 

-179 

-170 

124 

-112 

-97 

-89 

-84 

-83 

-70 

-5  0 

-28 

-15 

-7 

-5 

-2 

0 

12 

15 

15 

15 

38 

66 

73 

80 

81 

92 

111 

111 

132 

150 

159 

189 

201 

210 

235 

266 

285 

295 

321 

325 

351 

373 

394 

429 

446 

543 

549 

566 

552 

810 

£61 

1098 

MEAN  VALUE  =  80.812 
SDFV  =  410.026 
VIN  VALUE  =  -603 
MAX  VALUE  =  1971 
RANGE  =  2574 

1971 
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DATA  EILFNAME  PERSISTl 


■1005 
-454 
-248 
-214 
-143 
-110 
-89 
-46 

-26 

-17 

—  5 

-2 

33 

54 

51 

71 

93 

99 

159 

159 

207 

544 


-782 

-356 

-243 

-153 

-114 

-107 

-75 

-45 

-30 

-25 

-12 

-4 

20 

40 

56 

52 

61 

93 

115 

156 

177 

262 

555 


-577 

-257 

-214 

—  1 55 

-111 

-102 

-57 

-42 

-29 

-23 

-9 

30 
44 
59 
57 
92 
97 

153 
159 
190 
312 
1160 


MEAN"  VALUE  =  -4.565 
SEEV  =  279.752 
MIN  VALUE  =  -1005 
MAX  VALUE  =  1160 
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TATA  FILENAME  MODEL 2 


-579 
-53? 
-383 
-234 

-211 

-179 

-112 

-84 

-50 

-7 

0 

16 

61 
111 
159 

•21? 
255 
325 
394 
543 
652 

1098 


-560 

-422 

-272 

-234 

-211 

-170 

-97 

-S3 

-28 

-5 

12 

15 

73 

92 

132 

189 
235 
295 
351 
429 
549 
810 


•546 

■401 

■252 

■214 

■194 

•124 

-69 

-70 

-15 

-2 

15 

36 

60 

111 

150 

201 

266 

321 

376 

^46 

556 

c61 


MEAN  VALUE  =  52.60 
SDEV  =  335.425 
VIN  VALUE  =  -579 
MAX  VALUE  =  1098 
RANGE  =  1677 


61 


DATA    FILENAME    PERSIST2 


752 

-577 

41? 

-356 

248 

-243 

214 

-153 

143 

-114 

110 

-107 

-89 

-75 

-45 

-42 

-29 

-26 

-23 

-17 

-9 

-5 

-3 

-2 

30 

32 

4? 

44 

55 

59 

62 

57 

SI 

92 

93 

97 

115 

153 

j^CC 

169 

177 

190 

252 

312 

665 

MEAN 

VALUE    =   -11 

SDEV 

=   217.242 

MIN 

VALUE    =    -782 

MAX 

VALUE    =555 

RANG 

E    =    1447 

881 


-484 

-257 

-214 

-155 

-111 

-102 

-57 

-30 

-25 

-12 

-4 

20 

33 

54 

51 

71 

93 

99 

159 

159 

207 

544 
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EATA    FILENAME   M0DEL3 

*->  K*l   v_ 

-401 

-170 

-124 

-84 

-83 
-50 

15 

73 
150 

MEAN*    VALUE    =    -121.1 
SDEV    =   156.0  91 
MIN   VALUE    =    -538 
MAX    VALUE    =    150 
RANGE    =683 


TATA    FILENAME    PEESIST3 

-4S4 
-214 
-155 

-75 

-30 

-12 
-4 
62 

169 

169 

MEAN    VALUE    =    -57.6 
SDEV    =   183.947 
MIN   VALUE   =    -484 
MAX   VALUE    =169 
RANGE    =   653 
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TATA    FILENAME   n]0EEL4 

-211 

-179 

-69 

-5 
0 

12 

15 

38 

60 

SI 

92 
111 
201 
321 
325 
429 
652 
810 


MEAN  VALUE  =  149 
SDFV  =  261.680 
MIN  VALUE  =  -211 
VAX  VALUE  =  810 
RANSE  =  1021 


055 


54 


DATA  FILENAME  PERSIST4 

-243 
-114 

-29 

-25 

-23 

32 
40 

44 

54 

51 

6? 

71 

99 
159 
282 
544 
655 

MEAN  VALUE  =  93.278 
Sr-EV  =  2  29.622 
MIN  VALUE  =  -243 
M£X  VALUE  =  665 
PANGE  =  90S 
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TATA  FILENAME  M0EEL5 

-234 

-214 

-211 

-70 

-28 
-2 

210 
285 
295 
446 
549 
881 

MEAN  VALUE  =  158.917 
SEFV  =  333.733 
*IN  VALUE  =  -234 
MAX  VALUE  =881 
RANGE  =1115 


DATA  FILENAME  PERSISTS 

-248 

-110 

-102 

-57 

-45 

93 

93 

97 

115 

153 

190 

312 

MEAN  VALUE  =  40.917 
SDEV  =  148.878 
VIN  VALUE  =  -248 
MAX  VALUE  =  312 
RANGE  =  560 
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DATA  FILENAME  M0DEL6 

-579 
-562! 

-54:3 

-422 
-333 

-272 
-252 
-234 

-194 

-112 

-97 

-15 

-7 

16 

55 
111 
132 
159 
1S9 
235 
255 
351 
375 

— '  —  »_ 

555 
1098 

MEAN  VALUE  =  30.704 
SDEV  =  353.155 
MIN  VALUE  =  -579 
MAX  VALUE  =  109S 
RANGE  =  1577 
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.TA  E 

ILENAME  PIKS 

-762 

-577 

-419 

-356 

-257 

-214 

-163 

-143 

-111 

-127 

-59 

-42 

-25 

-17 

-9 

-3 

-2 

20 

30 

33 

56 

59 

61 

92 

169 

177 

207 

ISTo 

MEAN 

VALUE  =  -88. 

.519 

SLEV 

=  223.996 

MIN 

VALUE  =  -782 

MAX 

VALUE  =  207 

PANGE  =  969 

- 
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DATA  FILENAME  MODEL7 

-550 

-401 

-763 

-252 

-234 

-170 

-124 

16 

73 

60 

81 

92 

152 

159 

199 

210 

255 

325 

351 

394 

445 

543 


MEAN  VALUE  =  55.8^4 
SDEV  =  287.282 
*IN  VALUE  =  -550 
MAX  VALUE  =  543 
RANGE  =  1103 


DATA  FILENAME  PERSIST? 


-782 

-214 

-214 

-153 

—1  c5 

-114 

-111 

-110 

-76 

-5? 

-30 

-23 

-17 

-3 

20 

54 

62 

159 

169 

169 

177 

190 

MF/N 

VALUE  =  -48.591 

STEV 

=  203.037 

M I N 

VALUE  =  -782 

MAX 

VALUE  =  190 

RANGE  =  972 
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EAT  A    FILENAME    M0EEL8 

-579 
-548 

—  t^TQ 

-422 

-272 

-234 

-214 

-211 

-211 

-179 

-97 

-69 

-64 

-5  3 

-25 

-15 

-7 

-5 

-2 

12 

15 

15 

33 

56 

111 

201 

235 

295 

378 

549 

556 

652 

810 

881 

1098 


MEAN    VALUE    =    61. 
SDEV    =   390.058 
MIN    VALUE    =    -579 
KAX   VALUE    =    1098 
RANGE    =    1677 


086 
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DATA  FILENAME  PERSISTS 

-57? 

-484 

-418 

-355 

-257 

-243 

-107 

-102 

-89 

-45 

-42 

-29 

-25 

-12 

-9 

-5 

-4 

-2 

30 

32 

33 

44 

58 

59 

51 

81 

93 

93 

97 

115 

207 

282 

213 

544 

565 

MEAN  VALUE  =0.0 
SDEV  =  241.952 
MIN  VALUE  =  -577 
MAX  VALUE  =  555 
RANSE  =  1242 
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DATA 


FILENAME 

VCDEL9 

-194 

-112 

-83 

-70 

01 
■a 

111 

132 

285 

321 

429 

MEAN    VALUE    =    81.1 
SDEV    =    193. 22c 
fIN    VALUE    =    -194 
VAX    VALUE    =    429 
RANGF    =    623 


DATA 


FILENAME    P 

EESIST9 

-243 

-143 

-26 

42 

67 

71 

'92 

99 

153 

166 

MEAN 

VALUE 

=    27.3 

SDEV 

=    124. 

982 

MIN 

VALUE    = 

-246 

MAX 

VALUF    = 

166 

RANG 

E    =   416 

73 


DATA  FILENAME  MOLEL10 

-503 
-543 
-538 
-422 
-383 
-272 
-214 
-211 
-194 
-179 
-112 

-97 

-64 

-83 

-70 

-2 

0 

61 

•   159 

325 

579 

552 

810 

661 
1096 

MEAN  VALUE  =  22.92 
SDEV  =  450.575 
"KIN  VALUE  =  -60  3 
MAX  VALUE  =  1096 
F.ANGE  =  1701 
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DATA  FILENAME  PERSIST10 

-1005 
-57? 
-4.04 

-41S 
-257 
-24.3 
-143 
-111 

-99 

-26 

-23 

-12 
-2 
20 

51 

93 

93 

97 
153 
159 
168 
292 
544 
555 

MEAN  VALUE  =  -38.72 
SDEV  =  335.521 
MIN  VALUE  =  -1005 
MAX  VALUE  =  665 
RANGE  =  1670 
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DATA  FILINGS  M0CEL11 


-560 
-252 
-234 
-170 
_cg 

-29 

-7 

12 

15 

38 

73 

92 

111 

150 

201 

235 

285 

321 

378 

429 

543 

556 


-401 
-234 


-211 

-124 

-52 

-15 

-5 

15 

15 

66 

80 

111 

132 

189 

210 

265 

295 

351 

394 

446 

549 

1971 

MFAN  VALUE  =  140.023 
STEV  =  372.330 
MIN  VALUE  =  -560 
MAX  VALUE  =  1971 
FANC-E  =  2531 
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TATA  FILFNAME  PEHSIST11 


.792 

■24? 
■214 

■155 

■lie 

■102 

-57 
-42 
-29 

-17 

-3 
32 

42 

54 

52 

92 

115 

159 

192 
312 


-355 

-214 

— 167 

-114 

-127 

-75 

-45 

-32 

-25 

-9 

-4 

32 

33 

44 

58 

61 

71 

99 

159 

177 

227 

1180 


!*FAN  VALUE  =  7.83  5 
SEEV  =  243.521 
MIN  VALUE  =  -7=2 
MAX  VALUE  =  1132 
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APPENDIX  C 
HISTOGRAMS  OP  PREDICTION  ERRORS 

The  following  pages  are  histograms  of  the  prediction 
errors.  These  plots  show  the  distribution  of  each  data 
set.  The  number  of  tins  were  chosen  such  that  there  were 
at  least  two  date  points  per  bin  whenever  possible. 
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APPENDIX  T 
DELTA  *i  -  DELTA  ?  VALUES 

The  following  pages  contain  the  data  sets  used  in  the 
paired  sanple  t  test.  These  values  were  obtained  by 
subtracting   froj  Delta  v  the  corresponding  Delta  ?  value. 
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r^IA    FILENAME    PA  I  FEE  1 

-4:93  -423                               -325 

-30  7  -27? 

-235  -225 

-165  -130 

-109  -95 

-85  -84 

-57  -54 

-46  -20 

-G  -5 

24  27 

31  33 

4  0  40 

46  52 

71  89 

104  109 

145  166 

194  222 

237  256 

270  320 

352  35S 

397  402 

555  573 

791  1017 

MEAN    VALUE    =    106.594 
SEEV    =    304.525 
MIN    VALUE    =    -493 
MAX    VALUE    =   1156 
FAN'GE    =    1549 


-251 

-223 

-124 

-95 

-7  2 

-51 

_  1  * 

X  T 

17 

28 

1  Q 

44 

64 

97 

14  0 

19  0 

222 

267 

349 

ry  r-*   *■» 

3oo 

784 

1156 
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TATA  FILENAME  PAIREE2 


-403 

-270 

-225 

-130 

-96 

-94 

-54 

-23 

-6 

27 

33 

43 

52 

66 

133 

155 

222 

256 

320 

356 

402 

673 

1317 

■325 

■251 

■223 

■124 

-95 

-72 

-51 

-14 

17 

25 

39 

44 

54 

97 

140 

163 

222 

267 

349 

366 

533 

764 


-33^ 
-235 

-155 

-109 

-66 

-67 


-Q 

24 

31 

43 

4P 

71 

104 

145 

194 

237 

7P? 

352 

397 

565 

791 


MEAN  VALUE  =  99.561 
S2EV  =  271.398 
WIN  VALUE  =  -403 
MAX  VALUE  =  1017 
RANGE  =  142  0 
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TATA  FILENAME  PAIRED! 

-307 
-223 
-124 
-109 

-95 
l7 

190 

237 

256 

257 

402 

*J  <_'  <*J 

764 

M2A.M  VALUE  =  139.546 
SDIV  =  305.344 
MIN  VALUE  =  -307 
MAX  VALUE  =  7C,4 
RANGE  =  1091 


134 


DATA 

FILENAME  PAIREE4 

-403 

-272 

-235 

-225 

-155 

-130 

-85 

-51 

-14 

-S 

27 

28 

33 

40 

52 

97 

222 

270 

320 

349 

352 

355 

397 

565 

673 

791 

1017 

1155 

MEAN 

VALUE  =  184.571 

STEV 

=  377.412 

MI  N 

VALUE  =  -403 

MAX 

VALUE  =  1156 

RANG 

E  =  1559 

105 


TATA  FILENAME  PAIREE5 

-493 

-307 

-27? 

-225 

-155 

-130 

-109 

-95 

-94 

-72 

-54 

-46 

-20 

-14 

— S 

-6 

17 

24 

27 

28 

33 

40 

52 

64 

10S 

140 

145 

130 

237 

320 

349 

673 

784 

1017 

1156 

ME/N  VALUE  =   94.171 
SDEV  =  340.108 
MIN  VALUE  =  -493 
MAX  VALUE  =1156 

RANGE  =  1549 
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ATA  FILENAME 

PAiREie 

325 

-23  5 

223 

-124 

109 

-95 

-94 

-45 

-20 

-14 

-p. 

1? 

24 

27 

2° 

31 

33 

38 

40 

40 

44 

46 

52 

71 

5S 

97 

140 

180 

194 

222 

222 

237 

2*5 

2^7 

32? 

343 

352 

3-58 

366 

33? 

533 

555 

673 

731 

MEAN    VALUE    -    132 
SEHV    =    231.451 
MIN    7ALrJE    =   -325- 
MAX    VALUE    =   791 
RANGE    =    1115 
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INITIAL  DISTRIBUTION  LIST 


No. 


1.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22.314 

2.  Library,  Code  0142 

Naval  Postgraduate  School 
Monterey,  CA  93943 

3.  Professor  Kenneth  L.  Davidson,  Code  53Ds 
Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey.  CA  93943 

4.  Professor  P.  J.  Larson,  Code  55La 
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Naval  Postgraduate  School 
Monterey,  CA  93943 

5.  Professor  John  M.  Pouldry,  Code  ?3Po 
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Naval  Postgraduate  School 
Monterey,  CA  93943 

6.  Joint  EW  Center/SE 
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